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Abstract 

We consider Kirchhoff equations for vibrating bodies in any dimension in 
presence of a time-periodic external forcing with period 2ir/uj and ampli- 
tude e, both for Dirichlet and for space-periodic boundary conditions. 

We prove existence, regularity and local uniqueness of time-periodic 
solutions of period 2n/ui and order e, by means of a Nash-Moser iteration 
scheme. The results hold for parameters (o>,e) in Cantor sets having 
measure asymptotically full for s — » 0. | 

1 Introduction 

We consider the Kirchhoff equation 

(1) u tt - Au(l + J \Vu\ 2 dxJ = eg(x,t) x G Q, t G R 

where g is a time-periodic external forcing with period 2it/uj and amplitude e, 
and the displacement u : fi x K — > R is the unknown. We consider both Dirichlet 
boundary conditions 

(2) u(x,t) = Vx G dn, t G R 

where f2 C R d is a bounded, connected open set with smooth boundary, d > 1, 
and periodic boundary conditions on R d 

(3) t) = it(x + 2ttto, i) Vm e Z d , iel i , t G R 

where = (0, 27r) rf . 

Equation JT]) is a quasi-linear integro-partial differential equation having the 
structure of an infinite-dimensional Hamiltonian system, with time-depending 
Hamiltonian 

H(u,v)= I — cfo + / J dx + I dx) - J egu dx. 

It describes nonlinear forced vibrations of a d-dimensional body (in particular, 
a string for d — 1 and a membrane for d = 2). 
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This model has been proposed first in 1876 by Kirchhoff [23J in dimension 
one, without forcing terms, with Dirichlct boundary conditions, namely 

(4) u tt - u xx (l + J ul dx)=0, u(0,t) = u(ir,t) = 

to describe transversal free vibrations of a clamped string in which the depen- 
dence of the tension on the deformation cannot be neglected. Independently, 
Carrier [14j and Narasimha [31] rediscovered the same equation as a nonlinear 
approximation of the exact model for the stretched string. 

Kirchhoff equations have been studied by many authors from the point of 
view of the Cauchy problem 

u(x, 0) = uq(x) , Ut[x, 0) = U\{x) 

starting from the pioneering paper of Bernstein [B] . Both local and global exis- 
tence has been investigated, for initial data having Sobolev or analytic regularity. 
See for example [H [33J M\ U EUl HI 123 and the rich surveys [37] . 

In spite of the wide study for the Cauchy problem, to the best of our knowl- 
edge nothing is known about the existence of periodic solutions of Kirchhoff 
equations, except for the normal modes. 

Kirchhoff himself observed that equation ([4]), thanks to its special symmetry, 
possesses a sequence of normal modes, that is solutions of the form u(x,t) = 
Uj(t) sinjx, j — 1,2,... where Uj(t) is periodic. In general, normal modes are 
solutions of the form Uj (t)ipj (x) where ipj (x) is an eigenfunction of the Laplacian 
on fi. 

In presence of a forcing term g{x, t) this symmetry is broken and normal 
modes do not survive (except in the one-mode case g{x,t) = gj(t)(fj(x)). In- 
deed, decomposing u{x,t) = Uj(t)ipj(x) shows that all components uj(t) 
are coupled in the integral term \X7u\ 2 dx, and problem JT]) is equivalent to a 
system of infinitely many nonlinear coupled ODEs, namely 

u'frt) + Xfa® (l + ]T \lu 2 k (t)) = egj(t), j = 1,2,... 

k 

where g(x,t) = gj(t)ipj(x) and A| are the eigenvalues of the Laplacian on Q. 

In this paper we prove the existence of periodic solutions of ([T]) . We consider 
the amplitude e and the frequency to of the forcing term g as parameters of the 
problem. We prove that there exist periodic solutions of order e and period 2tt/uj 
when s is small and (e, uS) belong to a Cantor set which has positive measure, 
asymptotically full for e — > 0. We prove regularity estimates for the solutions, 
both in Sobolev and in analytic classes, and local uniqueness (see Theorem H] 
Remark Q] and Theorem [2] in Section [2]) . 

There are two main difficulties in looking for periodic solutions of |T]). The 
first one is the so-called "small divisors problem" , caused by resonances between 
the forcing frequency u> with its overtones and the eigenvalues A| of the Lapla- 
cian on CI. Such a problem arises in the inversion of the d'Alembert operator 
dtt — A, whose spectrum 

{-lu 2 1 2 + X 2 3 : ZeN, j = L2,...} 
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accumulates to zero for almost every u>. For this reason, (<9 tt — A) -1 cannot 
map, in general, a functional space in itself, but only in a larger space of less 
regular functions. This makes impossible the application of the standard implicit 
function theorem. 

The other difficulty is the presence of derivatives in the nonlinearity. In 
general, little is known about periodic solutions of equations of the form 

u tt - u xx = ef{x, i, u, u x ,u t ,u XXl u xt , u tt )- 

This problem has been studied by Rabinowitz [35] in presence of a dissipative 
term au t , a ^ 0, and frequency u> = 1; by Craig [J5] for pseudodifferential 
operators 

u t t - u xx = a(x)u + b(x, Icy 3 ?/) = , (3 < 1 ; 
by Bourgain [T3] in cases like Utt — u xx + pu + u\ — and, for quasiperiodic 

1/2 

solutions, [JO] utt — u xx = a(x)u + edj (h(x,u)). We remark that, in general, 
the presence of derivatives in the nonlinearity makes uncertain the existence of 
global (even not periodic) solutions, see for example the non-existence results 
in |21[ I24| for the equation utt — cl(u x )u xx = when a > 0, a(v) = 0(v p ) near 
0,p> 1. 

Our proof overcomes these two difficulties by means of a modified Newton's 
method in scales of Banach spaces, that is a Nash-Moser method. At each step 
of the Newton's iteration we impose some "non- resonance conditions" on the pa- 
rameter uj to control small divisors. For these non-resonance frequencies we can 
invert the linearised operator, which is a perturbation of the d'Alembertian, los- 
ing some amount of regularity. In this way we construct inductively a sequence 
of approximate solutions. The loss of regularity, which occurs at each step of 
the iteration, is overcome thanks to smoothing operators and to the high speed 
of convergence of the quadratic scheme. 

The application of Nash-Moser methods to infinite-dimensional dynamical 
systems having small divisors problems has been introduced in the Nineties 
by Craig, Wayne and Bourgain, in analytic or Gevray classes [TB I fl2 | [T3] , ITS]. 
Further developments are for example in [jJl El [32] • 

This technique, combined with Lyapunov- Schmidt reductions, is a flexible 
alternative with respect to KAM procedures 38, 25[[31]. In particular, currently 
available KAM methods seem not to apply to the quasi-linear problem UJ. 

Since we deal with not only analytic, but also finite order regularity, the 
scheme we use here differs from that in [3 [S] and it does not rely on analyticity 
assumptions. Such a procedure goes back directly to ideas of the original meth- 
ods of 29, 30, 39 and it is developed in [7]. Recently [8] this technique has made 
possible to prove the existence of periodic solutions of nonlinear wave equations 
for nonlinearities having only C k differentiability. We point out that some of 
the difficulties of [5] are not present here, thanks to the special symmetry of the 
Kirchhoff nonlinearity. Moreover, the roles played here by space and time are 
inverted with respect to [8] (U [5] . 

We remark that small divisors problems become more difficult in higher 
dimension. For this reason, not many works deal with such problems when the 
dimension is larger than one (e.g. [IU[U]). ^ n that case, indeed, Aj have a sub- 
linear growth, see ©. In general this causes further difficulties in the inversion 
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of the linearised operators. In the present case, however, the structure of the 
Kirchhoff nonlinear integral term makes possible the inversion in any dimension 
(Section E| . 

Finally, we note that in case of periodic boundary conditions (j3J zero is an 
eigenvalue of the Laplacian. As a consequence, we have to solve a space-average 
equation which is not present in the Dirichlet case (0), see (|57p in Section [71 



2 Functional setting and main results 

Let 2tt/lo be the period of g. We look for solutions u with the same period. 
Normalising the time t — > tot and rescaling u — > e 1 / 3 ^, |T]) becomes 

(5) lo 2 uu — Am = (i^Au J |Vu| 2 dx + g(x, t) 
where fi := e 2 / 3 and g,u are 27r-periodic. 

2.1 Case of Dirichlet boundary conditions 

Assume that <9f2 is C°°. Let A 2 , <fj(x), j — 1,2,... be the eigenvalues and 
eigenfunctions of the boundary-value problem 

{— Aipj — A 2 <fj in il 
ifj =0 on dfl 

with J Q ip 2 dx — 1 and Ai < A2 < .... Weyl's formula for the asymptotic 
distribution of the eigenvalues gives \j = 0(j 1 ' d ) as j — > 00, thus 

(6) Cj^ d < X 3 < C'j 1 ^ Vj = l,2,... 

for some positive C, C depending on the dimension d and on the domain f2 (see 
e.g. [Ml Vol. IV, XIII.15]). 

By expansion in the basis {(fj(x)}, we define the spaces 

K, s (f2) := {v(x) =Y i V jVi {x) : £ \v,\ 2 \f e 2 ^ < 00} 
j 3 

for s > 0, a > 0. Spaces Vb, s with er = are used in [4j. They are the 
domains of the fractional powers A s / 2 of the Laplace operator. See [H [5D] 
for a characterisation. For instance, Vb,2 = H 2 (Vl) n -ffg(f2). We note that if 
it 6 7 , a (n) then A k u G ^(ft) for all < fc < (s - l)/2. 

Spaces with s = are used in [3], where it is proved that U CT >oV^ ! o is the 
class of the (— A)-analytic functions, that is, by definition, the set of functions 
v(x) G Hq(SI) such that 



A k v G tfo(fi) and / vA fe vcfa; 



!/ 2 

<CA fe fc! Vfc = 0,l,. 



for some constants C, A In [3J it is observed that an important subset of 
U (T >oV r (T .o consists of the functions v{x), analytic on some neighbourhood of f2, 
such that 

A k v = on dtt Vfc = 0, 1, . . . 
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This subset coincides with the whole class of (— A)-analytic functions when dfl 
is a real analytic manifold of dimension (d — 1), leaving £1 on one side [27] . or 
when Q is a parallelepiped pQ. 

Clearly V a>a = {u G V afi ■ & s/2 u G V afi } and V aS) C V^, s C V a ^ for all 
s > 0, <7 > cr' > 0. Moreover, all finite sums J2j<N v j ( Pj( x ) belong to V a ^ s for 
all cr, s. 

We set the problem in the spaces X CTjS = iJ 1 (T, T4- )S ) of 27r-periodic functions 
u : T — > V^s, i i— > u(-,t) with ff 1 regularity, T := R/27rZ, namely 



:= {u(x,t) =J2 u j( t ) ( Pj( x ) ■ u i 6 ^(T.R), 

NIL :=EIK-Il^\ 2s e 2 ^<oo}. 



3>1 

Theorem 1. (Case of Dirichlet boundary conditions). Suppose that g G X ffiSo 
/or some cr > 0, So > 2d. Let s% G (1 + d, 1 + sq/2). There exist positive 
constants 5, C with the following properties. 

For every 7 G (0, Ai) there exists a Cantor set A 1 C (0, +00) x (0,07) 
of parameters such that for every {01,11) G Ay there exists a classical solution 
u(lu,/i) G X atSl of ([5]) @ . Sitc/i a solution satisfies 

IK w >/^)lksi < -C, ||w(a;,At)tt|| CT , sl -2 < -;C 
7 7^ 

and ji is unique in the ball {||u|| CT , Sl < 1}. 

The set A^ satisfies the following Lebesgue measure property: for every < 
uj\ < UJ2 < 00 there exists a constant C independent on 7 such that in the 
rectangular region TZ 7 :— (0)1,0)2) X (0,07) there holds 

'V^ 1 >i-a 7 . 

We recall that ([5]) is obtained from ([1]) by the normalisation t — > wr, and the 
rescaling u — > e 1 / 3 ^. Hence, going back, the solution u(o>, /1) of ([5]) found in 
Theorem [1] gives a solution of (fT]) of order e and period 2ir/oi. 

Remark 1. Theorem [T] covers both Sobolev and analytic cases: 

• (Sobolev regularity). If g belongs to the Sobolev space Xq^ Sq , then the 
solution u found in the theorem belongs to the Sobolev space ^o,si- 

• (Analytic regularity) . If g belongs to the analytic space X aa ^, then g G 
X aitSo for all o"i G (0, cr ). Indeed, 



Cso 

< 



exp[(cr - 0i)£] v (o"o - o-i)e 



V cr n - cti e / 



therefore 



3 

Since g G X ai _ So , the solution w found in the theorem belongs to the 
analytic space X aiiS1 C X ffli o- 
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Remark 2. If g{x, •) £ H r (T), r > 1, then the solution m of (JTJ) found in the 
theorem satisfies u(x, ■) € H r+2 (T) by bootstrap. 

Remark 3. (Nonplanar vibrations). We can consider the Kirchhoff equation 
for a string in the 3-dimensional space 

(7) u tt ~u xx (l + j \u x \ 2 dx^j = eg(x,t), g =(^ 9 ^j 1 u = (^J 

where the forcing g and the displacement u are ]R 2 -vectors belonging to the 
plane orthogonal to the rest position of the string, see [TH \5T\. In this case 
nonplanar vibrations of the string are permitted. 

Setting := + li M 2j||#i in the definition of the spaces X„ iS , 

Theorem [T] holds true for problem (J7J) as well. 

2.2 Case of periodic boundary conditions 

The eigenvalues and eigenfunctions of the Laplacian on T d are \m\ 2 , e m ' x for 
m G Z d . We consider a bijective numbering {mj : j G N} of Z d such that 
\rrij \ < |mj+i| for all j £ N = {0, 1, .. .}, and we denote 

X 2 := \ mj \ 2 , (p 3 {x) := e lm >- x Vj E N. 

We note that Ao = 0, (po{x) = 1 and Xj > 1 for all j > 1. Weyl's estimate © 
holds true for Xj as well, because the number of integer vectors m £ Z d such 
that |m| < A is 0{X d ) for A — > +oo (see [Ml Vol. IV, XIII.15]). 
We define 

X a , s := [u(x,t) =J2 u i( t ) { Pj(x) ■ Uj £ iT x (T, M), 

i>o 

Mis := ||«o|||ri +E INI^f ^ < °°}- 

Theorem 2. (Case of periodic boundary conditions). Suppose that g £ ^cr,s n 
/or some a > 0, s$ > 2d, and 

(8) / g(x,t)dxdt = 0. 

Let Si £ (1 + d, 1 + so/2). There exist positive constants 5, C with the following 
properties. 

For every 7 £ (0, 1) there exists a Cantor set A-y C (0, +00) x (0, 07) of 
parameters such that for every (uj, /i) G A 7 there exists a classical solution 
u(uj,fi) e X aySl of (J5]) (|3]) satisfying 

/ u(u, /i)(x,t) dxdt = . 

J(0,2TT) d + 1 

Such a solution satisfies 

(9) \\u(u,n)\\ a , ai <^(l + ^)c, |K^) tt || CT , Sl _ 2 < -^C 

7 V U) z / 7W 2 

and if is unique in the ball {J^ 27r y+i u (x, t) dxdt = 0, ||u|| . >Sl < 1}. 
The set A 1 satisfies the same measure property of Theorem^ 
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Remark 4. If u(u>, (i) is a solution of ©©, then also u(uj, fi) + c, c G K, solves 
©©. 

2.3 Outline of the proof 

The rest of the paper is devoted to the proof of the theorems. In Sections l3|4|5|6l 
we develop the details for the proof of Theorem [TJ then the same calculations 
are used to prove Theorem [5] in Section [7J 

In Section [3] we perform the Nash-Moser iteration to construct the approxi- 
mating sequence (u n ), for /i small and (w,/i) belonging to smaller and smaller 
"non-resonant" sets A n . Avoiding resonances allows to invert the linearised 
operator at each step of the iteration. 

In Section [J] we prove that u n converges to a solution of the Kirchhoff equa- 
tion if (w, /i) G A n for all n. Local uniqueness of the solution is also proved. 

In SectionOwe prove that the intersection of all A n is a nonempty set, which 
is very large in a Lebesgue measure sense. 

In Section [5] we prove the invertibility of the linearised operator for G 
A n and we give an estimate on the inverse operator. 

In Section [7J we complete the proof of Theorem [T] and we prove Theorem [5] 

3 The iteration scheme 

We fix a > once for all. In the following, we write in short X s := X a , s , 
\\u\\ s := UmIIo-^. We remark that all the following calculations holds true both in 
the Sobolev case a = and in the analytic case a > 0. Indeed, the only index 
used in the present Nash-Moser method is s. 

We set the iterative scheme in the Banach spaces X s endowed with the 
smoothing operators P„, defined in the following way. We consider a constant 
X G (1, 2) and denote 

(10) N n := exp( X ") 

for all n G N. We define the finite-dimensional space 

Xj<N n 

and indicate P n the projector onto X^ n > (truncation operator). For all s,a > 
there holds the smoothing properties 

(11) \\P n u\\ s+a < K\\u\\ s V«6l s 

(12) ||(J-P n )u|| a < N~ a \\u\\ s+a VueX s+a 

where / is the identity map. We denote 

L w :=^ 2 a tt -A, /(it) := Au J n \Vu\ 2 dx, 
F(u) := L^u - ^/(u) - ng 

so that © can be written as 

(13) F(u) = 0. 
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Note that / is not a composition operator, because of the presence of the in- 
tegral. The map / is cubic: indeed f(u) = A[u, u, u] where A is the three- 
linear map ^4[u,i>,u>] = Au Vw o Vro dx. Moreover, since the integral term 
J n \\7u\ 2 dx depends only on time, there holds 

f(u)eX^ ViiGl W . 
The quadratic remainder of / at u is 

(14) Q(u,h) := f(u + h)-f(u)-f(u)[h] 

= Au I \Vh\ 2 dx + Ah I (2WuoVh+\Wh\ 2 )dx. 
Jn Jn 

We observe that, if a(t) depends only on time, then 

\\a(t)u(x,t)\\ s < ||a|| ff i||«|| a 

(we omit a factor given by the algebra constant of i? 1 (T)). As a consequence, 
by Holder inequality it is easy to estimate ||/(u)|| s , ||/'(u)[/i]|| s and \\Q(u, h)\\ s . 

We adapt the Newton's scheme with smoothing operators P n to the special 
structure of problem ([5|). We will construct a sequence (u n ) defining 

(15) u := 0, u n+ i := u n - F' '(un)" 1 ^^ - nf{u n ) ~ nP n+ ig] 
provided the linearised operator 

F'{u n ) : h i-> L u h - nf(u n )[h] 

admits a bounded inverse F'iun)- 1 onl(" +1 >. In this inversion problem a small 
divisors difficulty arises. We will prove (Lemma Q} that F'(u n ) can be inverted 
if the parameters (uj,fi) belong to some "nonresonant" set A n+ i, defined as 
follows. First, 

A := (0,+oo) x (0,1). 
By induction, suppose we have constructed A n and u n . We denote 

(16) a n (t) := [ \Vu n \ 2 dx 7 

Jn 

we consider the Hill's eigenvalue problem 

{y"+p 2 (l + Lia n (t))y = 
\y(t)=y(t + 27r) 

and indicate (p| n ^) 2 its eigenvalues, I G N. For 

r>d, 7 g(0,Ai) 

we define 

(17) A n+1 := {(cj,[x) G A n : |u;^ n) - Xj\ > ^ VA, < N n+1 , I G n| . 
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Remark 5. Note that for all /x,ro the set A n (p) := {ui : £ A n } is open. 

Indeed, for every < u>i < C02 < 00 the intersection (wi,^) H A n (fi) is defined 
by means of finitely many strict inequalities (see ^ and (|51 [) ). 

We fix a positive constant R such that, if u £ X\ and ||u||i < R, then 
a(t) := /„ |Vu| 2 dx satisfies ||a||fl-i < f and Halloo < 1/2. 

Lemma 1. (Inversion of the linearised operator). There exist two universal 
constants K\, K[ with the following property. Let u £ X^ n ' with \\u\\i < R. Let 
(w,ju) e A n+1 . If 



Lemma 2. (Construction of the approximating sequence). Let g £ Y Sa , so > 
2d. Let t £ (0?, so/2). There exist a choice for x in the definition (jTUJ) and posi- 
tive constants K, &, do, with 6(2 — x) > r + 1, satisfying the following properties. 



(Induction step). Suppose we have constructed u\, ...,«„ by (|15|) for (tu,fJ-) £ 
A n , n > 1, where each A k +\ is defined by means of u k by (117|) . and uj, £ 
Suppose that fi/~/ < Sq. Let 

hk+i ■= u k+ i - u k . 
Suppose that for all k = 1, . . . , n there holds 



Proof. (First step). Since uq = and (us, /i) £ A\, by Lemma Q] F'(0) is 
invertible and (fT5j) defines 



(18) ^\\u\\l +x <K[, 

then F'(u) is invertible, F'(u)- 1 : Jf( n+1 ) -> and 

(19) ||*»~^llo < — \\h\\ T -i V^l^. 



Proof. In Section [HI 



□ 





(23) 



^(O)- 1 ^) + M (/ - P 1 )g] = iiLz'Pig . 



By (|19p . the inequality (f5U|) holds true provided 



(24) 



ifilML < Kexp(-bx), t-1<s . 
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(Induction step). To define h n+ i by (|15p . we have to verify the hypotheses 
of Lemma □ By (JTTJ and IJ3TJI 



(25) \\hk\\r+i <Nl +1 \\h k \\ a <KZ exp[(-6 + r + l) X " 

7 

because hk £ l' 1 '. Then 

+00 



(26) |K|| r+1 < £ ||^|| T+1 < ^ exp[(-6 + r + l)x k ] 

fe=i ' fe=i 

which is finite for 6 > r + 1 . Thus condition Q18p is verified provided 



(27) 6>t + 1, X 2 (^j C 2 <^, Q, := 5>xp[(-6 + r + l) X fe ] . 

' k>l 

Since ||u„||i < ||it n || T +i, by {26]) we have ||u n ||i < P provided 

(28) K-C <R. 

1 

Since (u>, //) G A n _|_i, we can apply Lemma [T] and we define ft. n +i according 
to the scheme, namely 

(29) h n+1 := -F'(u n )~ 1 [F(u n ) + fi{I - P n+ i)g] . 
By (Unj) we have 

(30) IK+1II0 < — \\F(u n ) + fi(I - P„+i).g|| r _i . 

7 



By construction (|T5|) , u n satisfies 

F'(u n -l)h n = -LuUn-i + ^/(tt n _i) + /iP„#. 

By Taylor expansion F(u„) = F(u n -i) +F'(u n -i)h n - fJ,Q(u n -i,h n ), where Q 
is defined in (Hill. Thus 



(31) F(Un) = -/*[(/ - P«)<? + OK-l,y 

and (PU|) gives 

(32) ||ftn+l||o < — \\(Pn+l -Pn)9 + Q{u n -X,h n )\\ T -i. 

7 

Now (P n+1 - P n ) 9 = (J — P„)P„+i<?, then by dTJ 

||(P„ +1 - P„)flllr-1 < "4 ||Pn + 1.9|!r-l+/3 < "4 ||<?|L 

for 

t-1 + P<s q , p>o. 
To estimate ||Q(u n _i, /i ra )|| T _i, we note that 

\\Au n ^ f |V/i n | 2 || T _! < ||u n _i|| T+ i||/i n ||? < K^C \\h n \\l 

JQ 7 
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by (f2"5)) . and \\h n \\i < AT n ||/i n || by (fTTj) . For the second term, recalling that 
2<t + 1, 

\\Ahn / V(2u n _i + /l n ) O Vh n \\r-1 < ||/l„|| r+ l||/l„||o||2u„_i + /l„|| 2 

< 2K^C \\h n \\ r+1 \\h n \\ 
7 

by (123, and \\h n \\ T+1 < N^ +1 \\h n \\ by (XT). Then 

IIQK-i, M||r-i < 3^ - CoK +1 \\h n \\ 2 . 

7 



As a consequence, (|22j) holds true provided 

(33) ^ii||9lk < ^exp(~6 X " +1 ) 
and 

(34) 3X X £ C N^WhX < \ exp(-6 X " +1 ) 

7 2 



Condition (|33[) is satisfied for 

2#i||<?k 



(35) /3 > b X , K> 



exp[(/3 - 6x)x] 
and, by (f2Tj) . condition (f3"4f is satisfied for 



(36) t(2 _ x)>T+1 , ;< { HiM } ^, 

Since 2r < sq, we can fix % £ (1, 2) so close to 1 that 

T-l + (r + l)^- < s - 

2 - X 

Now we fix b such that 

6(2 - x) > r + 1 , r - 1 + b X < s 

and then we fix (3 as 

P = S - T + 1 . 

So (|35|) and (|24|) are satisfied for K big enough, and we fix K in such a way. 
Then (j27|) .(|28 | and (|36|) are satisfied for /i/7 small enough. □ 

4 The solution 

Lemma 3. (Existence of a solution). Assume the hypotheses of Lemma\^ and 
suppose that (w,e) 6 A n /or all b£N. TTien i/ie sequence (u n ) constructed in 
Lemma\t^ converges in X T+ ± to :— X)fe>i hk- u oo is a solution of (|13p and 

(37) ||«oo||t+i < -c 
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for some C. Moreover, (u n )u converges to (itoo)tt in X T —\, 

(38) IKUooJttllr-l^-^C 

70^ 

so Woo is a classical solution of (J5|). 

Proof. By ((2*5j) . the series ^ fe ||/ifc|| r+ i converges, u„ converges to Uoo in A T+ i 
and (|37l) holds true. 

By (|23J) there holds uP{hi) tt = A/ii + ^P g. By ([TTJ> and (|2l|) 



if u 

||A/ii|| T _i+ A t||i 3 off||r-i < — exp[(-o + T + l)x]+/i||<?IU, 

7 

so that 

(39) ||(/ll)tt||r-l<C-^ 

for some C (recall that 7 < Ai). For n > 1, by (J29|) and (J3TJ) 

(40) F'(u n )/i n+ i = (i[(P n+ i - P n )g + Q(u n -i,h n )], 
thus 

Lo 2 {h n+1 ) tt = Ah n+1 + fi(f (u n )[h n+1 ] + (P„ - P n _x)g + Q(u n -x,h„)). 
By (ESI),®,® and (03} we get 

(41) \\(K +1 )u\\r-i <C-^exp[(-fo + T + l) X n+1 ]. 

It follows that (u n )tt converges in X r _i, (uoo)tt S Z T _i, so that has regu- 
larity H 3 C C 2 in time and ([55]) holds true. As a consequence F(u n ) converges 
to F(uoo) in X r _i. 

On the other hand, by ([31]) , J33]) and (|34ll 



II^K)llr-i<^exp(-6 X n+1 ), 

and P(u„) -> in X T _ X . Then F( Uoo ) = 0. □ 

Remark 6. We will prove in Lemma [5] that the set {(lj,/j) G A„ Vn € N} is 
nonempty and has positive, large measure. As a consequence, the sequence (u n ) 
of Lemma [3] is defined for all (w,/i) in that large set. 

Lemma 4. (Uniqueness of the solution). Assume the hypotheses of Lemma\^ 
There exists Si £ (0,o"o] such that, for /1/7 < 5\, Uoo is the unique solution of 
(Unj) in the ball {v € X T+1 : |M| T +i < !}■ 



Proof. Suppose v is another solution of (fT3[) . with ||w|| T+ i < 1. Let v n := P n v. 
Projecting the equation F(v) = on iW gives 

L^Un = ^t(/(u„) +P„5), R n (v) := Av n \V(v - v n )\ 2 dx. 
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Since u n solves (fl"5j) . that is 

F'(u n )h n+ i = -L^Un + (J,(f(u n ) + P n +ig), 
the difference w n := v n — u n satisfies 

LojWn - fJ,(f(v n ) - /(«„)) - ^R n (v) + //(P n +i - P n )g = F'(u n )h n+1 . 
Since f(v n ) - f(u n ) = f'(u n )[w n ] + Q{u n ,w n ), applying P'(w„) _1 

(42) w n = h n+1 +/iF'(ii„)" 1 [g(M„,w„) +Rn(v) - (P n+1 -P n )g] ■ 
Now, by (dH) 

\\liF'(u n )~ 1 Q(u n , w n )\\ < Ki - \\Q( 

7 

By assumption and (|3"T|) 

||^n||2 < ||W„|| T +1 < ||«]]t+1 + ||Uoo||r + l < C 

and by (g5]) ||u n || 2 < IKJIr+i < C, so that u>„)|| T _i < C||w„|| and 

\\^F'(u n y 1 Q(u n , w n )\\ <K\ — C ||tu n ||o < ^ ll^nllo 

7 -s 

provided /i/7 is small enough. Thus (|4"2"| gives 

\ IKIlo < ll^+l||o + #1 - IKW - (P„+l - P„)$||r-1 ■ 

2 7 

By (HH) and (fT2^) the right-hand side tends to as n — > 00, so that ||w n — tt n ||o — > 
0. Since v„ converges to v and u„ to Uoo in Xq, it follows that v = Uoo- D 

5 The Cantor set of parameters 

Lemma 5. (Regular dependence on the parameter to). Assume the hypotheses 
of Lemma\j^ There exist 82 6 (0, S{\ such that all the maps 

h n : A n n {(co, fi) : /x/7 < 5 2 } -> X (n) , (w, /x) i-> /i„(w, /x) 

are differentiable w.r.t. uj and 

(43) ||a««n||o<C-^ 
for some C . 

Proof. U\ = h\ 6 is defined for (w,/i) 6 Ai and it solves oj 2 (h\)tt — 

Ahi+fiPog. Recalling Remark[5]and Lemma[TJ by the classical implicit function 
theorem it follows that hi is differentiable w.r.t. u>. Differentiating w.r.t. lu gives 

2w(hi)tt + L u [d u hi] =0. 
We apply L" 1 and by (H|) and (131?]) 

||0J»l||o<C-£- 



13 



for some C. 

Assume that for n > 1 

(44) R,Mo<C-j- exp[(-6 + r + l)x fe ] Vfc = 1, 



/i„+i solves (|40|). then it is differentiable w.r.t. w. Differentiating (|40|) gives 

F' (u n )[d u h n+ i] = -2u(h n+1 ) tt + pbf" (u n )[d u u n , h n+1 ] + /j,d u (Q(u n -i, h n )). 
We apply F'^Un)^ 1 and observe that 

||F>„)-V(/>n+l)tt]||o < — IKAn+OttHr-l < exp[(-6 + r + 1) X " +1 ] 
by (gT]). To estimate 

(45) K 1 ^\\f"(u n )[d ul u n ,h n+1 } 
7 

we write all the integral terms and apply (TTTj) and (f5|) to each of them. We 
write here the calculations for two terms, the other ones are analogous. First, 



(46) \\A(d u u n ) / Vu„ o V/i n+ i|| T _i < \\d u u 

n || t+1 \ \ u n || 2 || ||o 

JO 

and ||9a,n n || T +i < ||e?u,u n ||o 

Rn„|| < E H^-Mo ^ C -T- E ex P[(- 6 + r + = HT 1 > 



IK|| 2 < |K||r+i < C^/7 by (25]), 

^ +1 H^+illo < ^ ex P [(-& + r + l) X n+1 ] by 

7 

so that (gBJ) < C(^/7) 2 (/i/7 2 w) exp[(-6 + r + As second example, 

||A/l n+ i / Vu„ o V(<9 u U n )|| T _l < ||/l„+l||o iV^+JllWnlUII^Unllo- 



'n+1 

ju 

At the end we have 



M <c(^Y-%- exp[(-6 + r+l)x" +1 ] 



j 3 

-7/ 7 'a; 

for some C, so that P%|) holds true for k = n + 1 provided /x/7 is small enough, 
independently on n. 

Finally, ||flL«n||o < ELi ll^fcllo and dm implies (031). □ 

Lemma 6. (The Cantor set). There exist S3 < 82 such that the Cantor set 
A 7 := Aoo l~l {(w, n) : fi < (S37}, 7 G (0, Ai), has the following measure property. 

For every interval I = (0)1,0)2) with < u>i < C02 < 00 there is a constant C 
depending on I such that, denoted by 1Z~/ the rectangular region TZ~ — lx(0, 63"/), 

l 1~L^ ("I «/4'-y 
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Proof. Wc fix /i, we recall that A n (/i) := {ui : (ui, fi) G A n } and define 

E n := A n (p) \ A n+1 {n), n e N. 

We have to prove that U ne jqE n has small measure. As a consequence, its com- 
plementary set Aoodi) := n ne nA n (n) will be a large set. Let 



We note that f2™ = for all j, n because 7 < Ai and — for I ~ 0. Suppose 
that u> £ Then CXj < ui < CXj for some C,C by JHI]). Moreover 

Q>n || 00 

by dill). Recalling the definition (TTBJ) of a n , by (03]) and 



|9 w o n || ff i = ||2 / Vw„o V(9 w u„)dx|| ff i < 2||m„|| 2 ||||i9 w u„||o < C ^r- 
Jn 7 w 



for some C, so that 



By (JSTJ it follows that 



l^ (n) |<c;4- 



(47) 



» 



/i 3 i 

^- > 7 

7 ,3 a; 4 



provided /i/7 is small enough, say /i/7 < ^3. Fix < u>% < L02 < 00. If 
f2™ ; (~l (0)1,0)2) is nonempty, then 



(48) m< ^ <Caaj ^, ^(-A,,-^.)^) 
for some C, C", C" . Since £Jq = Ua^a^;^^;, we have the estimate 

(49) " " 



\E Q n(Q u u 2 )\< ]T E l^l<7C* £ ^ 



Aj<iVi ieA(j) 

for some (7 depending on (Qi,Q 2 ). 

To estimate \E n n (<Di,(Zi 2 )|, n > 1, we notice that 



A 3 '<AT„ +1 , 1>1 

For the sets with N n < Xj < N n+ i we use (14*51) and we get 

N„<Aj<N„ +1 JV„<Aj<JV„ + i ? 
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where C is the constant of To estimate the remaining sets, suppose that 

oj € ny, for some Xj < N n , I > 1. Then by ([55j) 

|A,-«P l (n_1) | < lA.-^l+cb^-^l 
T 

< — + 2w/^||a„ - a„_i||oo . 
A i 

Since < CA^ and \\a„ - a n -i\\ H i < ||/i n ||o||2u n _i + h n \\ 2 , by (JH]) , (T35]) we 
have 

|A 3 ■ -^"^l < J exp(-6 X n ). 

i ' 

Thus 

,3 



and by (|i7|) 



|ny,,rii4 n 0i)| <Cwa4 ex p(- & x")- 
It follows that 

ieA(j) 

for some C depending on (u>i, ui 2 ). Now, by (O Aj < 2V W implies j < (N n /C) d , 
then 

r (N n /C) d 

E A J ^ £ - C " / - 

Xj<N n j<{N,JC) d 

for some C" . As a consequence 

E ^exp(-b X n )<Cexp[(-b + d+l)x n ]. 

Aj<Af„ 



Then 



|B„n(a>i,« 2 )|<7C( E ^+exp[(-b + d+l) X n }) 



N n <\ j <N„ +1 3 

for some C , and 



U £„n(^ 2 )| < 7 c(E^ + E ex p[(- fe+d+1 )^] 



The first series converges because by © 
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being t > d. The second series converges because b > d + 1. Thus 

| (J s„n (0)1,0)2) | <c 7 

and the relative measure of Aoo(fi) in (0)1,0)2) satisfies 

l^oo (m) n (0)1,0)2)1 > x _ c 
0)2—0)1 

for some C. 

Finally, we integrate in in the interval where fi/j < 63, 
\A 00 nRy\= \Aooifi) n (0)1,0)2)1 d/x. 



□ 

6 Inversion of the linearised operator 

In this section we prove Lemma H Let u <E l'"U £ The linearised 

operator is 

F'(u)h — L u h — fj,f'(u)[h] =Dh + Sh 
where we split F'(u) in a diagonal part 

Dh := uj 2 h tt - Ah(l + (i J |Vu| 2 da;) 

and a "projection" part 

Sh:=—[iAu / 2Vuo V/irfs. 
Jn 

We recall here some results on Hill's problems. The proof is in the Appendix. 

Lemma 7. (Hill's problems). Let a(t) be 2% -periodic and ||a||oo < 1/2- The 
eigenvalues p 2 of the periodic problem 



(50) 

form a sequence {p 2 }i S N such that 



y" +p 2 (l + a(t))y = 
y(t) = y(t + 2tt) 



(51) \l<Pi<2l VI e N. 

For a, (3 2tt -periodic, \\a\\oo, \\P\\oo < I/ 2 , 

(52) \p l (a)-pt(p)\<2l\\a-l3\\ O0 VI 6 N. 

TTie eigenfunctions ipi(t) of (|50p /orm an orthonormal basis of L 2 (T) w.r.t. 
the scalar product 

(u, w) £ 2 = / m;(l + a) dt 



17 



and also an orthogonal basis of H l (T) w.r.t. the scalar product 



2 K 



(u,v) H i = / u'v' dt + (u, v) L 2 . 
J a 

The corresponding norms are equivalent to the standard Sobolev norms, 

(53) ~\\y\\ L * < \\y\\ Ll <2||y|| £ », \ \\y\\ H , < \\y\\ K < 2 

Lemma 8. (Inversion of D). Let u e jW, a(<) := J Q |Vu| 2 da; £ iT^T), 
< 1, ||a||oo < 1/2. Let p 2 be the eigenvalues of the Hill's problem 



(54) 



y" +p 2 (l + pa(t))y = Q 
y{t) = y(t + 2tt) . 



// (lo, fi) satisfy the non-resonant condition 

\upi - Xj | > VAj- < N n+ i, I EE, 

then D is invertible, D^ 1 : X^+V -> X^ n +^ and 

(55) P _1 %< -\\h\W-i V/ ie x(" +1 ) 

7 

/or some constant C . 

Proof. If h — Y^, hjit) <fij(x), then Dh — Y Djhj(t) fj(x), where 
Djz(t) = u?z"(t) + X 2 z(t)p(t), p(t) := 1 + pa(t). 
Using the eigenfunctions ?/>;(£) of (|54l as a basis of i? 1 (T), 

D^(t)=53(A|-wV)^^(*)p(*)» * = £3iV'i(*). 

ZGN ieN 

and Kj := (l/p)Dj is the diagonal operator {X 2 — uj 2 p 2 }i e -^. Since \X 2 —uj 2 pf\ > 
1 /X T f 1 for all Xj < N n+1 , we have that Kj is invertible and 



ZGN J 



By (ED H^-^Hhi < 4(AJ~V7) 11*11^- Since DT^ = Kj\zfp) and 
is smaller than a universal constant, 

ry 

Since L> _1 /?, = £\ Dj 1 hj(t)^j(x) we obtain (55]). □ 

Lemma 9. (Control of S). For all s > 0, if u E X s+ 2 then S : Xq — > X s is 
bounded and 

||5/i|| s </i|H| s+2 ||u]] 2 ||/i||o VhEX . 
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Proof. Since f n V« o Vh dx does not depend on x, 



Ait / Vu o V/i dx 


< \\Au\\ s 


/ Vu o Vhdx 




Jn 


s 


Jn 





J^VuoVhdx = Y.jXjUjitfhjit), so H/nVuoV/idajllirifD < ||u|| a ||fe||o by 
Holder inequality. □ 

Proof of Lemma [Q F'(u) = D + S = (I + SD- r )D where I is the identity 
map. Since satisfies (|55p . we have to prove the invertibility of I + SD~ X 
in norm || || T -i- By Neumann series it is sufficient to show that 



(56) 

By Lemmas [5] and H] 



SD^hWr-i < i \\h\\ T -i Vft g X (n+1) . 



\SD-'h\\ T ^ x < m||u|| t _ 1+2 ||u|| 2 ||D-^|| < ||u||^ +1 ||/i|| T _i 



because 1 1 ] ] 2 < ||w|| r +i. Thus the condition 



implies ^ and by Neumann series + SD -1 ) -1 /i|| T _i < 2||/i|| T _i. 



□ 



7 Proof of the theorems 

Proof of Theorem [TJ Let g 6 X G:So and 2d < 2(s\ — 1) < sq as assumed in 
the theorem. We apply Lemma [2] with 

t := s% — 1 . 

The construction of the sequence (u n ) is possible provided the parameters (w, (i) 
belong to A n for all n G N. Lemma H] assures that, for /x/7 sufficiently small, 
the set Ay of parameters satisfying this property is a nonempty set, which is 
very large in a Lebesgue measure sense. Lemmas [3] and |4] complete the proof. 

□ 

Proof of Theorem [2l Let g G ^cr,s with 2d < 2(si — 1) < so as assumed 
in the theorem. We consider a Lyapunov-Schmidt reduction splitting the space 
X a:S in two subspaces X a ^ s — Y (W D ^cr, s ), 

F:={y(t) Gtf^R)}, W := {« G l ,o : «>(a:, i) = ^ ic^t) ^(x)} . 

We denote Hy, H\v the projectors on Y, W, and observe that Hy is the map 

u 1 ► / u(x, i) dx. 

J(0,27r) d 

We define 

g (t) := Hy£/, 5(2;, t) := n^S- 
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We decompose u(x, t) = y(t) + w(x, t), y £Y, w £ W, and note that 

/(«) = f(y + w) = f(w) e w. 

Then projecting equation on Y gives 

(57) u> 2 y"(t) = ngo(t) (Y equation) 
while projecting it on W 

(58) L u w — fx(f(w) + g) (W equation). 

Equation ([57]) is an ODE. With direct calculations (or by Fourier series) we see 
that ((57)) admits 27r-periodic solutions if and only if 



(59) / g o (t)dt = 

Jo 

and (|59"|) is just assumption ©. We note that, if y(t) solves (j57|). then also 
y(t) + c solves ([57]) . for all c 6 1. Moreover, the unique solution y(t) of ([57]) 
such that J Q v y(t) dt = satisfies 

\\y\\m < \\v"\\m<±\\ go \\ Hl . 

To solve ([58]). we consider all the calculations in Sections I3|4|5I6I replacing 
X a:S with X a ^ s (~l W and Xj , <pj (x) with Xj , ipj (x) , j > I. It follows the existence 
of a unique solution w £ X a , si n W of ([58]) satisfying 

Then u = y + w solves ([5])(j3|). Since 

\\u\\l, si = \\yf Hl + \\w\\l, Sl , \\u tt \\l Sl - 2 = \\y"\\ 2 m + \\wu\\l, Sl - 2 , 
we obtain estimates ©. □ 

8 Appendix 

Proof of Lemma The proof follows from classical results in [TH [22]. First, 
if y" + p 2 (l + a)y = 0, then 

y n dt=p 2 (1 + a)y 2 dt , 
o Jo 

so that p 2 > because 1 + a is positive, p 2 , — is an eigenvalue, the cor- 
responding eigenfunctions are the constants, and all the other eigenvalues are 
positive. 

By [ini Theorem 2.2.2, p. 23], for every k £ N both P2k+i and p2fc+2 satisfy 
\(k + lf <p 2 < 2(fc + l) 2 
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and (|51[) follows. 

The equivalence (|53p and the orthogonality of (tpi) w.r.t. ( , )#i can be ver- 
ified by direct calculations. 
To prove (l5"2"l) . we define 



g (0)(i) = g (t) := 1 + o(i) + - a(t)), G [0, 1] 

and the "Liouville's change of variable" t — > £ 



1 /■* m . 1 



« = /(?) ^ £ = := - J V?(«) d«, c := 2^l V «(*) 



We note that p 2 , satisfy 
(60) | 
if and only if p 2 , z(£) satisfy 
(61) 
where 



y"(t)+p 2 q(t)y(t) =0 
t/(t)=y(t + 27r) 



/'(0 + c a [p 2 -Q(/(0)MO = 
z(£) = z(£ + 2tt) 



The operators T(i?) : z i-» — z" + c 2 Q(f)z are selfadjoint in L 2 (0, 2w). We apply 
[221 Theorem 3.9, VH.3.5, p. 392] to the holomorphic family {T(i?) : i? G [0, 1]} 
(see [12 Definition VII.2.1, p. 375, Example 2.12, VH.2.3, p. 380 and Example 
6.13, III. 6. 8, p. 187] to verify the hypotheses of the Theorem in the present case) 
to prove that the eigenvalues and eigenfunctions of (frTTj) are analytic in As 
a consequence, the eigenvalues and eigenfunctions of (|60|) are analytic in d as 
well. This allows us to differentiate the equation 

+Pi(0) 2 (l + a + - a))tM*) = 
w.r.t. 1?. Recalling that qipf dt — 1, multiplying by ^>j(i?) and integrating 

/-2tt 

2 J 

Since < 2? and q(i9) > 1/2, 



toWl < « / ' ^nrr^ sW;W 2 * < - "lie 



g(0) 

and ^ 

|pi(/9)-pl(a)|< / |^p z (7?)|dT?< 2/||/3 -a| 



□ 
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